CHARACTERISTIC POLYNOMIALS 
OF SAMPLE COVARIANCE MATRICES 



H. KOSTERS 

■ Abstract. We investigate the second-order correlation function of the charac- 

' teristic polynomial of a sample covariance matrix. Starting from an explicit 

formula for a generating function, we re-obtain several well-known kernels from 
random matrix theory. 



in 



1. Introduction 



Ph ' Characteristic polynomials of random matrices have found considerable attention 

. in recent years, one reason being that their correlations seem to reflect the corre- 

; lations of the eigenvalues [BHTI IBH2[ IFSI IBDSI EEl EE EHl EEl iKoTl lKo2] . 

In this article we investigate the second-order correlation function of the charac- 
teristic polynomial of a sample covariance matrix. 

, Complex Sample Covariance Matrices. Let Q be a distribution on the 

[ real line with expectation 0, variance 1/2 and finite fourth moment b, and for 

' given n,m £ N with n > m, let X := X{n,in) := {Xij)i=i^^^^^n;j=i,...,m denote 

. the n X m matrix whose entries Xij are i.i.d. complex random variables whose 

I real and imaginary parts are independent, each with distribution Q. Let X* = 

' X*(n, m) denote the conjugate transpose of Z. Then the Hermitian mxm matrix 

^ ■ Z := Z{n,m) := (l/n) X{n,m)* X{n,m) is called the (complex) sample covariance 
matrix associated with the distribution Q. For G C, let 

^ \ f{n, m; fi, u) := E {det{X{n, m)*X{n, m) — fi) det(X(n, m)*X(n, m) — u)) 

d denote the second-order correlation function of the characteristic polynomial of the 

"unrescaled" sample covariance matrix. 

Real Sample Covariance Matrices. Let Q be a distribution on the real 
line with expectation 0, variance 1 and finite fourth moment b, and for given 
n,m G N with n > m, let X := X{n,m) := (^ij)t=i,...,n;j=i,...,m denote the n x m 
matrix whose entries Xij are i.i.d. real random variables with distribution Q. Let 
= X''" {71,171) denote the transpose of X. Then the symmetric mxm matrix 
Z := Z{n,m) := (l/n) X{n,m)^ X{n,m) is called the (real) sample covariance 
matrix associated with the distribution Q. For /i, 1/ G C, let 

/(n, m; //, u) := E (det(X(n, m)'^X{n, m) — //) det(X(n, m)^X{n, m) — v)) 

denote the second-order correlation function of the characteristic polynomial of the 
"unrescaled" sample covariance matrix. 
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In the special case where Q is the Gaussian distribution, the random matrix 
Z is also called a (complex or real) Wishart matrix, since its distribution is the 
(complex or real) Wishart distribution (see Anderson [An] or Muirhead [Mu]). 

It will always be clear from the context whether we are considering the complex 
or real case. More precisely, we are interested in the asymptotic behavior of 
the values f{nj\f,m]\f;fi]\f,h'j\f) as N ^ oo, where njy := A^, tun := N — a for 
some fixed natural number a, and the parameters /i^v and vn are finally used to 
"zoom in" at certain interesting regions of the spectrum of the (unrescaled) sample 
covariance matrix, namely the bulk of the spectrum, the soft edge of the spectrum, 
and the hard edge of the spectrum. (See Section 3 for details.) 

In the complex setting, we will recover, in all the regions previously mentioned, 
the well-known kernels for the "correlation functions" of the eigenvalues, namely 
the sine kernel, the Airy kernel and the Bessel kernel. (See Section 3 for details.) 
Thus, our results indicate that these kernels are universal in that they arise in 
connection with quite general sample covariance matrices (as described at the 
beginning of this section), albeit at the level of the characteristic polynomial. 
It is conjectured that this universality also holds at the level of the eigenvalue 
themselves, but so far this conjecture has only been proven for a restricted class 
of sample covariance matrices in the bulk of the spectrum (see Ben Arous and 
Peche ^BPj ) and under stronger moment conditions than ours for the edges of 
the spectrum (see SOSHNIKOV [So] and Tao and Vu [TV], respectively). 

Similar results will also be obtained in the real setting. However, the results 
for the correlation functions of the characteristic polynomial are somewhat more 
different from the results for the eigenvalues here. 

To obtain our results, we follow the strategy proposed by Gotze and Kosters 
|GK] in the context of Wigner matrices. First, we obtain an explicit expression for 
an exponential- type generating function of the second-order correlation function 
of the characteristic polynomial. Second, we recover the well-known kernels from 
random matrix theory by asymptotic analysis. 



In this section we derive the generating function of the second-order correlation 
function of the characteristic polynomial. Since the derivation is very similar for 
complex and real sample covariance matrices, we present the details only for the 
complex setting and restrict ourselves to a few comments for the real setting. 

The following well-known representation will prove useful: 

Lemma 2.1. Let X e M{n x m,C). Then, for X^O, 



2. Generating Functions 




Proof. This is an immediate consequence of the matrix factorizations 




In 

CA-^ D - CA-^B 
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^ B \ _ ( ^- BD-^C BD-^ \ ( 
C D ) - [ Irn )[C D 

where A G M{n x n, C), B £ M{n x m, C), C G A4(m x n, C), G A4(m x m, C), 
and D and ^ are invertible. □ 

A matrix of the form , where X is a (complex) random matrix with i.i.d. 

entries, is also called a (complex) chiral matrix. Thus, Lemma 12.11 establishes 
a connection between the characteristic polynomial of a sample covariance matrix 
and that of a chiral matrix. Chiral matrices seem more convenient for our purposes, 
as they give rise to neater row and column expansions. 

To illustrate our approach, we start with the first moment of the characteristic 
polynomial, which has already been determined by Forrester and Gamburd 
|FGj by means of a combinatorial argument. Fix A G C, and put 

where X := X{n,m) is a matrix of size n x m with i.i.d. matrix entries satisfying 
our standing moment conditions. Note that this definition is meaningful even if 
n = or m = 0; we then have f{0,m) = A™ and /(n,0) = A". (In particular, 
the determinant of the 0x0 matrix is defined to be 1.) 

Let us derive a recursive equation for /(n, m). Suppose that m > 0. Then, doing 
a row and column expansion about the last row and the last column, we have 

det( ^ =Adet^ 



m\\* 




X* Mm J V (^'" 

+ X:(-ir^^-^An.^.,n.det( f;-"] ], (2.1) 

where an upper index [i : j] indicates that the ith row and the jth column of the 
corresponding matrix are deleted. Taking expectations, using independence and 
noting that 

(as follows from our standing moment conditions), we obtain 

/(n, m) = Xf{n, m — 1) — nf{n — 1, m — 1) . (2-2) 

Observe that this is meaningful even if n = 0, since the second term vanishes in 
this case. Hence, recalling the identity 

/(n,0) = A", (2.3) 



we see that the values /(n, m) are completely determined by recursion over m. 

nomials 



For arbitrary a G C, the Laguerre polynomials L^\x) are defined by 
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(see e.g. Equation (5.1.6) in SzEGO [SzJ). It is well-known that the Laguerre poly- 
nomials satisfy the relations 

mL(--"^Hx) = -xLtr'\x) + nLtT\x) , (2.4) 

Lt'^\x) = l, (2.5) 

(see e.g. Equations (5.1.10) and (5.1.14) in SzEGO [Sz]), through which they are 
completely determined by recursion over m. 

Comparing (|2.2p . (12. 3p with (|2.4p . (12. 5|) . it is easy to see (by induction on m) 
that for any n,m >0, 

f{n,m) = (-l)'"m!A"~"4r"^)(A2). (2.6) 

Hence, using Lemma |2.H we find that the first moment of the characteristic poly- 
nomial of a (complex) sample covariance matrix satisfies 

Kdet{XIm- X*X) = (-l)™m!L;;-'"(A). (2.7) 



In fact, this was already proved by Forrester and Gamburd |FG] by means of 
a combinatorial argument. 

Remark. Equation (12. 7p remains true for real sample covariance matrices (under 
the respective moment conditions). 

Remark. We have obtained a recursive equation over m for the moments /(n, m) 
Of course, for symmetry reasons, it is clear that it is equally possible to derive 
a recursive equation over n. Indeed, starting with a row and column expansion 
about the first row and the first column, we obtain 

/(n, m) = Xf{n — 1, m) — mf{n — 1, m — 1) . 

instead of (j2.2p . together with the initial condition 

/(0,m) = A-. 

Since the Laguerre polynomials satisfy the relation 

L\^'^\x) = {-xr/m\, 

(see e.g. Equations (5.1.13) and (5.2.1) in SzEGO [Sz]), they can also be calculated 
by recursion over n, and (12. 6p follows (by induction on n). 

Let us now turn to the second moment of the characteristic polynomial. Similarly 
as above, we first consider 
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( X* J„ ) ( X* ut 



/^det ( ^^fi]^, ^^J" ) + X] (-l)^+^-^X,,^X,,^det 



Note that due to our standing moment assumptions, we have 

EXij = 0, EXfj = 0, E|Xy|2 = l, E|Xijf = 26^ + 1 , 
Hence, expanding the product and taking expectations, we obtain 

/(n,m)=^.Edet^^^,^„,^, J ( (Xl ' H)* ,/_J 

- ,Y.E\X,,J det^ Mn-J'^*l (Xl-H)* J 

- .gE|X,.| det^ ^/„-J^^*l(^[-¥ J 
+ ^E|X,„| det ^ J det ^ J 
+ j;E|X,„p|X,,„pdet ( ) det 



To shorten notation, let us introduce the auxiliary functions 



f(n,m,a, ■):=E (det ( ' u 



det 



where a = (ai, a2; «3; 04) G {'jIj^}^ (a dot representing no deletion), 
as well as the abbreviations • :=(•,-,•,•), 01 :=(•,-, 1, 1), 10 := (1, 1, • , • ), 
11a := (1,1,2,2), 11b := (1,2,1,2), 11c := (1,2,2,1). (No other values of a 
will be needed. Moreover, the value 11^ will be needed for the real setting only.) 
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Note that f{n, m, ■ , ■) = f{n, m). We then have 



f{n,m) = nuf{n,m- 1) 

— iinf{n, m — 1, 01, • ) 

— vnf[n, m — 1, 10, • ) 

+ (26 + i)ra/(n- l,m- 1) 
+ n{n - l)f{n, m-l, 11a, ■ ) 

+ n(n- l)/(n,m- 1,11c, •)• (2-8) 

Note that the preceding formula was obtained by starting with an expansion 
"in the m-dimension" . Of course, it is possible to derive an analogous formula by 
starting with an expansion "in the n-dimension" ; we then have 

/(n, m) = iJiffin — 1, m) 

— fj,mf{n — l,m, • , 01) 

— h'raf[n — l,m, • , 10) 

+ (26+ l)mf{n - l,rn - 1) 
+ m(m — l)/(n — 1, m, • , 11a) 

+ m(m- l)/(n- l,m, -.llc), (2.9) 
where f{n,m, ■ ,a) is defined by 



/(n,m,.,a):=Eldet( ] ) (Xl'-l)* ./J 



0!4:o;3j 



and the possible values • , 01, 10, llyi, Hb, He for a are the same as before. 

Here is a list of recursive relations which can be derived using the above arguments: 

f{n, m, • , ■ ) = nuf{n, m - 1, • , • ) 

— finf{n, m — 1, 01, • ) 

— i'nf{n, m — 1, 10, • ) 

+ (26+i)n/(n-l,m-l, •, •) 
+ n(n- l)/(n,m- 1,11^, •) 

+ n(n- l)/(n,m- 1,11^, •) (n>0,m>l), (2.10) 

/(ra, m, • , 01) = {+n)f{n, m - 1, • , • ) 

-n/(n,m- 1,10, •) (n>0,m>l), (2.11) 

/(n, m, • , 10) = {+u)f{n, m - 1, • , • ) 

-ra/(n,m- 1,01, •) (n>0,m>l), (2.12) 



Characteristic Polynomials of Sample Covariance Matrices 7 

f{n, m, ■ , 11^) = nuf{n, m - 2, • , • ) 

— finf{n, m — 2, 01, • ) 

— i^nf{n, m — 2,10, • ) 
+ n/(n — 1, m — 2, • , • ) 

+ n(n- l)/(n,m- 2, 11^, •) (n>0,m>2), (2.13) 

/(n, m, • , 11*^) = n/(n - 1, m - 2, • , • ) 

+ n(n- l)/(n,m- 2, 11'^, •) (n>0,m>2), (2.14) 

/(n, m, • , • ) = MJ^/l'^ - 1, m-, • , • ) 

— /im/(n — l,m, ■ , 01) 

— umf{n — 1, m, • , 10) 

+ (26 + i)m/(n- l,m- 1, •, •) 
+ m(m — l)f{n — 1, m, • , 11"^) 

+ m(m-l)/(n-l,m, -.ll*^) (n>l,m>0), (2.15) 

/(n, m, 01, • ) = {+fJ')f{n - 1, m, • , • ) 

-m/(n- l,m, -,10) (n>l,m>0), (2.16) 

/(n, m, 10, • ) = {+v)f{n - 1, m, • , • ) 

— m/(n - l,m, • ,01) (ri>l,m>0), (2.17) 

/(n, m, 11^, • ) = /ui//(n - 2, m, • , • ) 

— jimf{n — 2, m, • , 01) 

— vmf{n — 2,m, • , 10) 
+ mf{n — 2, m — 1, • , • ) 

+ m(m- l)/(n- 2, m, -,11^) (n>2,m>0), (2.18) 

/(n, m, 11*^, • ) = mf{n - 2, m - 1, • , • ) 

+ m(m- l)/(n- 2, m, •,11'^) (n>2,m>0). (2.19) 

(When evaluating these recursive relations, there may arise some undefined terms 
with negative arguments, but this poses no problem as these terms are always 
accompanied by the factor zero. Similar remarks apply to the formulas below.) 
Together with the initial conditions 

fin, 0) = ifii^r , /(O, m) = if,ur , (2.20) 

the equations (I2.10p - (I2.19P determine the values f{n,m) completely. 

We will now derive a recursive equation involving the values f{n,m) only. 
By (|2.18p and (j2.19p with m replaced by m — 1, we have, for n > 2, m > 1, 

/(n,m-l,ll^, •) + /(n,m-l,ll^, •) 
= jiuf[n — 2, m — 1) — /i(m — l)f{n — 2,m — 1, ■ , 01) — iy{m — l)f{n — 2,m — 1, ■ , 10) 

+ 2(m - l)/(n - 2, m - 2) + (m - l)(m - 2) (^/(n - 2, m - 1, • , 11^) + /(n - 2, m - 1, • , 11*^)) . 
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By ()2.15p with n, m replaced by n — 1, m — 1, the two summands involving 11"^, 11 
satisfy, for n > 2, m > 1, 



(m- l)(m-2)(^/(n-2,m- 1, • , 11^) + /(n - 2, m - 1, •,ll'^)j 
= —fivfln — 2, m — 1) + fi{m — l)/(n — 2, m — 1, • , 01) + ^{m — l)f{n — 2, m — 1, 
+ /(n - 1, m - 1) - (26 + - l)/(n - 2, m - 2) , 

whence, for n > 2, m > 1, 

/(n,m-l,ll^, •)+/(^,m-l,ll^, •) 

= 2(m - l)/(n - 2, m - 2) + /(n - 1, m - 1) - (26 + i)(m - l)/(n - 2, m - 2) . 

(2.21) 

Plugging this into (12.10p and rearranging terms, we find that, for m > 1, 

f{n, m) = ^uf{n, m — 1) — finf{n, m — 1, 01, • ) — vnf{n, m — 1, 10, • ) 
+ (n(n + l)+n(26-|))/(n-l,m-l)-(n(n-l)(m-l)(26-|))/(n-2,m-2). 



(For n = and n = 1, p. 220 follows directly from p.lOp .) Thus, we have eliminated 
the terms involving 11^ and 11*-". 

The terms involving 01 and 10 can be eliminated by a similar substitution. 
Using ()2.16p and (I2.17|) with m replaced by m — 1, ()2.1ip and ()2.12p with n,m 
replaced by n — l,m — 1 and finally (I2.22p with n,m replaced by n — l,m — 1, 
we have, for n > 1, m > 1, 



= — (/U^ + v^)f{n — 1, m — 1) + (m — 1) ( /u/(n — 1, m — 1, • , 10) + vfin — 1, m — 1, 



= -{i? + v'^)f{n - 1, m - 1) + 2^lv{m - l)/(n - 1, m - 2) 

- (n - l)(m - 1) (/x/(n - 1, m - 2, 01, • ) + ^f{n - 1, m - 2, 10, • ) 



= — (//^ + v'^)f{n — 1, m — 1) + ijLv{m — l)/(n — 1, m — 2) 
+ (m — l)f{n — 1, m — 1) 

- (n(n- l)(m- 1) + (n - l)(m - 1)(26 - |))/(n - 2, m - 2) 
+ ((n - l)(n - 2)(m - l)(m - 2)(26 - |))/(n - 3,m - 3) . 

Plugging this into (j2.22p and rearranging terms, it follows that, for m > 1, 

/(n, m) = n(n + m)/(n — 1, m — 1) — n^(n — l)(m — l)f{n — 2, m — 2) 



+ iJLvf{n, m — 1) + ^vn{m — l)f{n — 1, m — 2) — (/i^ + v'^)nf{n — 1, m — 1) . 



(2.22) 




+ (26 - |)n(/(n - 1, m - 1) - 2(n - l)(m - l)/(n - 2, m - 2) 



+ (n - l)(n - 2)(m - l)(m - 2)/(n - 3, m - 3) 



(2.23) 
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(For n = 0, (|2.23p follows immediately from (j2.22p .) By symmetry, we also have, 
for n > 1, 

/(n, m) = (n + m)mf{n — 1, m — 1) — (n — l)m^(m — l)f{n — 2, m — 2) 
+ (26 - |)m(^/(n - 1, m - 1) - 2(n - l)(m - l)/(n - 2, m - 2) 

+ (n - l)(n - 2)(m - l)(m - 2)/(n - 3, m - 3)^ 

+ fiuf{n — 1, m) + )Ui/(n — l)m/(n — 2, m — 1) — (^^ + u'^)mf{n — 1, m — 1) . 

(2.24) 

Hence, the values f{n,m) may be computed recursively starting from (I2.20|) and 
using either <E1M or lEl^ . 

From now on, for n > 0, m > 0, let f{n,m) denote the second-order correlation 
function of the characteristic polynomial of the (complex) sample covariance matrix 
as defined in the Introduction. (For the rest of this section, we will usually omit 
the parameters ^, u, which are regarded as fixed.) Then, by Lemma l2.lt we have 
to make the replacements f{n,m) i-^ f{n,m) ■ (/uz/)""™, ^ ^jl, v"^ ^ v in 
the preceding two equations, thereby obtaining 

/(n, m) = n(n + m)/(n — 1, m — 1) — n'^{n — l)(m — l)/(n — 2, m — 2) 
+ (26- |)n(^/(n - l,m - 1) - 2(n - l)(m - l)/(n - 2,m - 2) 

+ (n - l)(n - 2)(m - l)(m - 2)/(n - 3, m - 3)^ 
+ /ii^/(n, m — 1) + fiun^m — l)/(n — 1, m — 2) — (/i + v)nf{n — 1, m — 1) 

(2.25) 

for m > 1 and 

/(n, m) = (n + m)mf{n — 1, m — 1) — (n — l)m?[m — l)f{n — 2, m — 2) 
+ (26 - |)m(^/(n - 1, m - 1) - 2(n - l)(m - l)/(n - 2, m - 2) 

+ (n - l)(n - 2)(m - l)(m - 2)/(n - 3, m - 3)^ 
+ /(n — 1, m) + (n — l)m/(n — 2, m — 1) — (/i + i')mf{n — 1, m — 1) 

(2.26) 

for n > 1, respectively, along with the initial conditions 

/(n,0) = l and /(O, m) = (/iz^)™ . (2.27) 

(Observe that in contrast to the equations for the chiral matrix, the equations for 
the sample covariance matrix are not fully symmetrical in n and m.) 

Let a := n — m denote the difference of n and m and assume that a > 0. 
(The case n < m could be reduced to this case by exchanging the roles of n and m 
and by multiplying with the appropriate power of fiv.) We will determine, for any 
a E N, the generating function 

^ f{m + am) ^ (2.28) 
^ (m + a)lml 

where \z\ < 1. 
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E{z) := exp ( -(^ + ly)- + (26 - |)z 



Put 

z 

We will show that for any a G N, the generating function (j2.28p is given by 

FM) ■■= E{z) (1 - TFtStIT^I - • ^2.29) 

fc=o ^' • 

To begin with, observe that Fa{z) defines an analytic function on the unit disc in 
the complex plane, with power series representation 

OO 

F^{z) = Y,c^{m)z^ , (2.30) 

»n=0 

say. Next, differentiating (|2.29p w.r.t. z, we obtain 



k=l 

Since 



1 2 I- z + 2z 1 + z 

+ 



z 1 — z z{l — z) z{l — z) 
it follows that 



(2.31) 

or (equivalently) 

F^(z) = (2zF;(z) + (2 + a)F„(z)) - {z^F'^^iz) + (2 + q)zF„(z)) 

+ (26 - |)(1 - z)2f«(z) - (^ + i/)F„(z) + //i/(l + z)Fo,+i{z) . 

(2.32) 

In terms of the power series coefficients Ca{m) defined by p.30p . this translates 
into the recursive relation, for m > 1, 

mca{m) = (2m + a)ca{m — 1) — (m + a)co(m — 2) 

+ (26 - I) {ca{m - 1) - 2ca(m - 2) + Ca(m - 3) 

+ fiuca+i{m - 1) + fivca+iim -2) - {fi + v)ca{m - 1) , (2.33) 

where terms with a negative argument are to be regarded as zero. For m = 0, 
we clearly have 

Cc,(0) = 1/q!. (2.34) 

Comparing ([2:25]) . ([2:271) and ([OS]) . ([23i]) . it is easy to see that the coefficients 
Ca(Tn) satisfy the same recursive equations as the values f{m+a, m)/((m+a)! m!) . 
Thus, as all values are uniquely determined by these recursive equations, it follows 
that the generating function (|2.28p is given by ()2.29p . for any q G N. 
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Let us summarize our result for the complex setting as follows: 

Proposition 2.2. For any a G N, the second-order correlation function f{n,m) 
of the characteristic polynomial of an unrescaled complex sample covariance matrix 
satisfies 

/("^ + «^"-^) _ p 

^ (m + a)!m! " ^ 



where b* := 2(6- |). 

Also, let us state the analogous result for the real setting: 

Proposition 2.3. For any a G N, The second-order correlation function f{n,m) 
of the characteristic polynomial of an unrescaled real sample covariance matrix 
satisfies 



^ (m + a)!m! " ^ 



m=0 



)oo 
I — n 



exp (-(;. + + 6*z ) ^^^^d - .)-^^— ^ , (2.36) 



where b* := b — 3. 



As already mentioned at the beginning of this section, the derivation in the real 
setting is essentially the same as that in the complex setting. That is why we do 
not give the full details of the proof of Proposition 12. 3|. but only mention some 
noteworthy changes: 

(i) In the real case, we have 

EXij = , EXfj = 1 , EXfj = b . 

In particular, as the second moment is not equal to zero anymore, we get an extra 
term n(n — l)f{n,m — 1, 11b, ■ ) in Equation (j2.8p . and an extra term m{m — 1) 
f{n — l,m, • , 11b) in Equation (|2.9|) . Also, the factor (26 + i) must be replaced 
with the factor 6 everywhere. 

(ii) Similar changes arise in Equations (|2.10p - (j2.19p . The recursive equation 
for the llB-terms is the same as that for the llc-terms, except that all occurrences 
of 11(7 must be replaced with occurrences of 11b- In feet, it is not hard to see 
that /(n, m, 11b, • ) = fin, m, 11c, ■ ) and /(n, m, • , 11b) = f{n, m, ■ , 11c) for ah 
n,m > 0. 

(iii) The analogue of (j2.2ip reads 

/(n,m- 1,11^, •)+/(n,m-l,ll^, . ) + /(n, m - 1, 11^, • ) 
= 3(m - l)/(n - 2, m - 2) + /(n - 1, m - 1) - 6(m - l)/(n - 2, m - 2) . 
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Consequently, the analogue of (|2.22p reads 

/(n, m) = /uz//(n, m — 1) — finf(n, m — 1, 01, • ) — vnf{n, m — 1, 10, • ) 
+ (n(n + 2) + n{h - 3))f{n - 1, m - 1) - (n(n - l)(m - l){b - 3)) /(n - 2, m - 2) , 

and the analogue of (j2.23p reads 

/(n, m) = n{n + m + l)f{n — 1, m — 1) — (n + l)n(n — l){m — l)/(n — 2, m — 2) 
+ (6 - 3)n(^/(n - 1, m - 1) - 2(n - l)(m - l)/(n - 2, m - 2) 

+ (n - l)(n - 2)(m - l)(m - 2)/(n - 3, m - 3)) 

+ /iJ^/(n, m — 1) + /ii/n(m — l)f{n — 1, m — 2) — (^u^ + u'^)nf(n — 1, m — 1) . 

Similar modifications are required for equations (j2.23p - (j2.26p . 

(iv) Defining Fa{z) as in Proposition ()2.3p and denoting by Ca{m) the coefficients 
in the power series representation Fa{z) = Ylm=o it follows by the same 

arguments as for (j2.33p that for m > 1, 

mca{m) = (2m + a + l)ca{m — 1) — (m + a + l)ca{m — 2) 

+ (6 - 3) (c„(m - 1) - 2c„(m - 2) + CQ(m - 3)) 

+ nuca+ii'm - 1) + fJ-J^Ca+iim -2) - {fi + u)ca{m - 1) . 

Similarly as in the complex setting, this is the same recursive relation as that for 
the values f{m + a, m) /{{m + a)\ m!) . 



3. Outline of the Asymptotic Analysis 

From now on, we will always assume that n = tin := N and m = rriN := N — a 
depend on N, where a is a fixed natural number. In particular, 171^/12^ < 1 and 
mAr/riTv — > 1 as ^ cxd. 

Let us begin with the complex setting. It is well-known from random matrix 
theory that under the above assumptions, the spectrum of the properly rescaled 
sample covariance matrix ^X*X is asymptotically concentrated on the interval 
[0,4], with distribution given by the Marcenko-Pastur density 

9(0 = ^7^(4^ (ee(0,4)). (3.1) 

Moreover, it is well-known from random matrix theory that the following 3 regions 
deserve particular attention: 

(i) the bulk of the spectrum, 

i.e. the region around a point ^ G (0,4), 

(ii) the soft edge of the spectrum, 
i.e. the region around the point ^ = 4, 

(iii) the hard edge of the spectrum, 
i.e. the region around the point = 0. 
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More precisely, it is widely expected that in these regions, the correlation function 
of the eigenvalues (see Mehta [Me] or Forrester [E]) is asymptotically given 
(after the appropriate rescaling so that the mean spacing between the eigenvalues 
is of order 1) by 



(i) the sine kernel 



(ii) the Airy kernel 



sin7r(x — y) , , 

§{x,y):=^ f^, 3.2 

7r(x - y) 



x-y 

(iii) the Bessel kernel 

^a{x,y):= ^- ^- ^ ^ — . (3.4) 

2(x - y) 

Typically, these results were first obtained for Wishart matrices (i.e., Gaussian 
sample covariance matrices) and then extended to more general sample covariance 
matrices: 

(i) For the bulk of the spectrum, Ben Arous and Peche |BPj estabhshed the 
emergence of the sine kernel for sample covariance matrices X*X such that 
the distributions of the entries of the matrix X are Gaussian convolutions. 

(ii) For the soft edge of the spectrum, Soshnikov [Soj established the emergence 
of the Airy kernel by means of the method of moments and a sophisticated 
comparison with the Gaussian case. 

(iii) For the hard edge of the spectrum, Tao and Vu [TVj established the emergence 
of the Bessel kernel (in the special case where the matrix X is a square matrix), 
also via a sophisticated reduction to the Gaussian case. 

It is the purpose of the following sections to show that the same kernels show up 
at the level of the characteristic polynomial, for quite general sample covariance 
matrices. Thus, our results add some support to the universality conjecture that 
the kernels always occur, irrespective of the choice of the underlying distribution Q. 
Also, our proofs are comparatively simple, and they are based on rather weak 
moment conditions, considerably weaker ones than those for the above-mentioned 
results from the literature. 

The starting point for our asymptotic analysis will also be the integral represen- 
tation 



f{n,m;n,u) ^ J_ 
n! ml 2-Ki 



j exp {-{^x + z^) + ^*^) 



oo 

which follows from (j2.35p by Cauchy's formula. Here, b* is defined as in (j2.35p . 
a := n — m, and 7 denotes a counterclockwise circular path of radius R < 1, 
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i? ~ 1 around the origin. (The precise choice of the radius will be specified in 
the following sections.) Since 



1/2 



^^^{k + ay.k\^' -y 1-z J 1^ 1 

where la denotes the modified Bessel function of order a, (|3.5|) may be rewritten 
as 

1 . exp(-(^ + .)^ + 



f{n,m;n,u) 



(3.6) 



n\ ml 2TTi (1 - z)2 . ((^i/^)i/2)« ^^+1 

We will always use the parametrization ^{t) := i?e**, — vr < t < +vr, and we will 
always assume that ^/z is chosen such that arg ^/z = on the positive real axis and 
arg-y/z is continuous otherwise, so argy^ G [— 7r/2, +7r/2] as z traverses the path 
7(t) := -Re**, — vr < t < +tt. Similar remarks apply to the choice of the branch of 
the modified Bessel function. 

In the course of our asymptotic analysis, we will use the following classical results 
about the asymptotic behavior of the modified Bessel function: For fixed a > 0, 

Uz) = ^^gS (1 + Oe{l/z)) {\z\ ^ oo, I argzl < f - e) (3.7) 



|/„(^)| <C^^Pfe^ (|z|^oo,|argz|<f) (3.8) 
\z\ 



and 



where C > is an absolute constant. (See e.g. Chapter 7 in Olver [01].) In fact, 
our results could even be extended to the more general situation where m]\j/ni\; 
tends to some constant 7 G (0, 1) as — > 00. However, this requires a more careful 
treatment of remainder terms (particularly the remainder terms originating from 
the modified Bessel functions) and will therefore not be pursued here. 

We will show that in all the above-mentioned cases, given the appropriate 
specialization of the shift parameters /x, v and the radius R, the main contribution 
to the integral comes from a small neighborhood of the point z = R and that 
we asymptotically end up with the above-mentioned kernels from random matrix 
theory. 

Furthermore, by essentially the same proofs, we obtain similar results for real 
sample covariance matrices. In this setting the correlation function of the charac- 
teristic polynomial is asymptotically given (after the appropriate rescaling) by 

(i) the "differentiated" sine kernel 

7r(x — [x — y)^ 

(ii) the "differentiated" Airy kernel 

2 Ai(x) Ai'(y) - 2 k\'{x) Ai(y) 



A(x,y) := 



{x - y) 



3 



, (x + j/)Ai(x)Ai(j/)-2Ai^(x)Ai^(j/) 
(x - yY 
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(hi) the "differentiated" Bessel kernel 



x,y) :-- 



4{x - yf 

These kernels are obtained from the corresponding kernels for the complex setting 
by applying the differential operators D := — (in cases (i) and (ii)) 

and D := {n Sy ~ ^ Sx) (^^ ^^^^ (iii))- (Besides that, to obtain the above 
formulas for the "differentiated" kernels, we have to use the differential equation 
for the Airy and Bessel function.) Since the derivation of the results for real 
sample covariance matrices is essentially the same as for complex sample covariance 
matrices, we will refrain from giving the details of the proofs, but only state the 
final results. 



4. ASYMPTOTICS IN THE BULK OF THE SPECTRUM 

In order to zoom in around a point in the bulk of the spectrum, we have to make 
the replacements 

fJ- ^ N£, + ^/g{£,) and v ^ N$,+v/g{$) 

for some ^ G (0,4), i/ S M. (Recall that our sample covariance matrices are 
unrescaled and that g{£^) denotes density of the Marcenko-Pastur distribution.) 
Moreover, we have to multiply by the proper scaling factor in order to obtain 
a non-degenerate limit. We then have the following result: 

Theorem 4.1. Let f denote the second-order correlation function of a complex 
sample covariance matrix satisfying our standing moment conditions. For any 
a G N, ^ G (0, 4), fi,u £R, setting 

Zn{^, fi, v) := (iVV + Ni{^l + //I/)"/' exp {-Ni - i(^ + v)) , 



we have 



hm [g{i)-'ZM{i,is/g{0,y/gm , , 

N~*cxi \ n\ ml 

= exp(6*)S(/i,z.), (4.1) 

where b* := 2(6 — |) and § is defined as in (j3.2p . 
Proof of Theorem \4-l\ By (|3.6p , we have 
f{n,m;NC + fi,NC + iy) 



n! ml 

1 r exp (-(iVe + Ul^ + i^))^, + b*z) . {w{z)) 
2^i X (1 - zf ■ z("+™)/2 T ' ^^'^^ 
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where we have put 

w{z) = wn{z) := 2y^N^e+N^{fi + u)+f,iy 

for abbreviation. We take R := 1 — We reduce the proof of Theorem 14.11 

to the following two claims: 

On the one hand, we will show that, for any (fixed) a > 0, 

1 r exp (-(iVe + Ul^ + i^))t^ + b*z) ■ Ia{w{z)) 

I ^ ' — > Sia) , (4.3) 



(1 _ zY ■ z("+'")/2 



where 7' = j'{a) denotes the restriction of the path 7 to the set {—a/N,+a/N), 
and 



S{a) :-- 



exp( 



(b*) exp((l - iO • (1 - iu) - lifi - u^-'/il - ^n)) 



4^3/2^1/2 



(1 - iM)3/2 



du . 



On the other hand, we will show that for any d > 0, there exists some ao((5) > 
such that for a > ao{5), we have 



1 

2TTi 



exp (-(iVe + ^(/i + i^))^ + b*z) ■ Ia{w{z)) 



(1 _ 2)2 . z{n+m)/2 



dz 
z 



< 5 



(4.4) 



for all G N large enough. Here, 7" = 7" (a) denotes the restriction of the path 7 
to the set (-vr, +7r) \ {-a/N, +a/N). 

It is easy to see that the main result (I4.ip follows by combining ()4.3p and (14. 4p 
and by using Laplace inversion. Indeed, firstly, (|4.3p and (14.4jl imply that 



1 . exp (-(iVe + ^(/^ + ^))^ + &*^) • ^(^(^)) dz 



N-*oo 



where 



S{oo) :-- 



(1 - Z)2 . ^(™+m)/2 2 

exp(6*) /■^+*~ exp((l - i^)^ - i(/^ - ^^)TV^) 



5(00) , 



47^3/2^^1/2 



1— ioo 



,3/2 



dz . 



Secondly, by Laplace inversion, for t > 0, a € 



1 
27ri 



l+ioo 



' 1— ioo 

(see e.g. p. 245 in [Er] ) . whence 

5(00) = iexp(6*) ((1-10/0 



,3/2 



1/2 



dz 



sm 



2 sin(a-v/t) 



1/2 



(/z - i.) ((1 - i^/e) 



1/2 



Replacing the local shift parameters n, v with iij giX)-, v I g^X)-, dividing by gi^) and 
noting that ((1 - = 7rg{0 yields (gl]). 



as 
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To prove (|4.3p and (j4.4p , it turns out convenient to rewrite the integrand in (j4.2p 

fl{^)f2{z) 



where 

and 
Since 

and 



(1 - Z)2 2{n+m)/2 ' 

fi{z) := exp (^-(iVe + ^(/^ + ^))r3^ + b*^) ■ (+"'(^)) 
f2iz) := exp (-w(z)) • la {w{z)) . 



w 



1-^ l-z~ (1 + Vi)(l- V^) ~ 1 + 

we find that 



/i(^) = exp ( - ive[^ - + ^)[^ - - ^)V(A^e) 



1 - z 
+ 6*z + C'(l/iV" 



1 -z 

Here and in the sequel, we take the convention that the imphcit constants in the 
O-terms may depend on a, h* fi, u (which are regarded as fixed). 



We first prove (j4.3p . In doing so, we use the notation Oa to denote a bound 
involving an implicit constant depending also on a (in addition to a, 6*, ^, i/). 
Substituting z = i?e** and t = u/N on the left-hand side in (j4.3p . we obtain 

1 h{l{u/N))h{l{u/N)) 
2vriV (1 - -i{u/N)Y ■ (7(n/iV)){"+'")/2 ' ^ ■ > 

Now, by straightforward Taylor expansion, we have the approximations, for \u\ < a, 
hiRe'"'") = exp ( - 1«1 - m) - 1(M - >^fr^l(l - in) + b') (1 + O.(JV-l)) 

(1 - i?e^"/^)2 = (1 - m)ViV^ (1 + Oa{N-^)) , 

^^gm/A^)(m+n)/2 ^ g^p ( - (1 - m)) (l + Oa{N~^)) . 

For the second approximation, we have used the observation that for |n| < a, 
u;(i?e*"/^) = 2N'^i/{l - iu) {l + Oa{N-^)) , as well as the asymptotic approxi- 
mation (j3.7p for the modified Bessel function. Inserting the preceding approxi- 
mations into (j4.5p . it follows that the left-hand side in (j4.5p converges to 

exp(6*) [+'^ exp((l - \0 . {I - iu) - - 1^)2^7(1 " iu)) 



which proves (j4.3p . 
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We now turn to the proof of (j4.4p . It suffices to show that for a > ao(5), 



1 

2^ 



[-TT,~a/N] U 



\h{Re^')\\f2{Re^')\ 



dt<5. 



(4.6) 



Write = re*"^ with r = ^Jl - 1/iV, 99 G [-7r/2, +7r/2]. Then, 



Re 



1 - re'^f 



> 



1 



1 + re*"^ / 1 + r2 + 2r cos (/3 AN 



1 - re*"^ 



1 + re*<^ 



re 



1 - (re*'^)2 



|1 -re*^| 
|1 + re*"^! 

1 

< 



< 1 



N, 



from which it follows that 



h{Re'') = 0{l) . 



(4.7) 



Also, for N G'N sufficiently large. 



\w{Re'^)\ > 



re 



up 



1 - {re^^y 



1 



> hN£ 1 TT > jNe 



and 



it\2 



>o. 



Re {w{Re^')) = 2^ N^e + + u) + ^iu ■ Re iyY^^-^ 

from which it follows using ()3.8p that 

|/2(i?e**)| = O (|1 -i?e*Y/V^/^) • 
Finally, 



(4.8) 



Thus, by ()4.7p . (14. Sp . (14. 9p and symmetry, it remains to show that for a > aQ{6), 



1 



1 



dt < 6N^/^ 



(4.10) 



27r |l-i?e^*|3/2 

To this end, let c > be a small constant such that cost < 1 — c^t'^ for all |t| < vr. 
Then we have 

|1 _ i?e**p = (1 - Rcostf + (i?sint)2 = (1 - Rf + 2i?(l - cost) > c^t^ 
and therefore 

/ \l - Re'Y^/^ dt < {ct)-^/^ dt = 2ica/N)-^/\ (4.11) 



which entails ()4.10p by choosing aQ{6) > large enough. 
This completes the proof of Theorem 14. li 



□ 
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With essentially the same proof, we obtain the following result: 

Theorem 4.2. Let f denote the second-order correlation function of a real sample 
covariance matrix satisfying our standing moment conditions. For any a G N, 
^ G (0,4), E M, setting 

Zn{^, fi, ly) := {N^f + N({fi + + iivf''' exp (-iV^ - + z.)) , 
we have 

N-^oo V n! m! / 

= exp(6*)S(/x,i/), (4.12) 

where b* := {b — 3) and § is defined as in i\3.9\} . 



5. ASYMPTOTICS AT THE SOFT EDGE OF THE SPECTRUM 

In order to zoom in at the soft edge of the spectrum, we have to make the 
replacements 

lj.^m + and v ^ AN + 

where ^u, i/ G M. (The factor 2^/^ is for convenience.) We then have the following 
result: 

Theorem 5.1. Let f denote the second-order correlation function of a complex 
sample covariance matrix satisfying our standing moment conditions. For any 
a G N, ^, 1/ G M, setting 

Z^ifi, u) := (^16N^ + 4(/u + u)N^/^ + fiuN^/^Y^^ exp (^-4iV + i^)iV^/^) , 

we have 

N^oo \ n'.ml I 

= exp(6*) A(^,zy) , (5.1) 

where b* := 2(6 — |) and A is defined as in ()3.3p . 
Proof of Theorem \5.1[ By (|3.6p , we have 

1/3 /(n,m;4iV + //iVV3^4jV + i/jVV3) 

n! m! 

^1/3 I- exp (-(4iV + ^(/x + z.)iVi/3) 1+1 + 5*^^ . 
~ 27rz 
where now 



X (1 - Z)2 . ^(n+m)/2 T ' ^^'^^ 



t(;(z) := 2\ 16N'^ + 4(^ + u)^/^ + /izyiV2/3 



1 - z 

This time we take R := 1 — Similarly as for the bulk of the spectrum, 

we reduce the proof to the following two claims: 
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On the one hand, we will show that for any (fixed) a > 0, 
Ari/3 r exp (-(47V + i(^ + i.)iVi/3) i+£ + 5*^^ . /„(«;(^)) 

(5.3) 

where 7' = j'{a) denotes the restriction of the path 7 to the set {—a/N^/^,+a/N^^^), 
and 

exp(6*) /■+" exp (^(1 - iuf - + - iu) - ^{fi - - iu)-') ^ 

""^^^■■=^^1^ (Y^W 

On the other hand, we will show that for any 5 > 0, there exists some ao{6) > 
such that for a > ao((5), we have 



< S (5.4) 



^1/3 . exp (-(4iV + i(/z + z.)iVi/3) i+£ + . /^(^„(^)) 

27ri 7 ,/ (1 - z)2 . 2{n+m)/2 J 



for all G N large enough. Here 7" = 7" (a) denotes the restriction of the path 7 
to the set (-^,+7r) \ {-a/N^/^,+a/N'^/''^). 

Using the integral representation for the Airy kernel given in |Ko2j , the theorem 
may then be deduced using similar arguments as in the previous section. Indeed, 
first of all, ([53]) and dEU imply that 

^1/3 r exp (-{AN + i(^ + z.)iVV3) i±| + t*z) ■ Uw{z)) 

' ^ > A{oo) , 



27ri (1 - z)2 . ^(n+m)/2 

where 

Aioc) = "^P(^*) exp(^.3_i(^^^)^_^(^_^)2^-i) 



Substituting z = 2^/^^, dz = 2'^^^dz, and shifting the path back to the line 1 + 
(which is easily justified by Cauchy's theorem), we obtain 

2^/3 exp (6*) /-W"" exp {^z^ - |(^ + ^^) 2^/^z - - uf 



16^3/2^ J^_^^ ^3/2 

Making the replacements ji ^ 2'^/^^^ v 1— > 2^/3 a,nd multiplying by 2^/3^ 
we further obtain 

exp(6*) /-W"" exp (^^3 _ ^ ^) _ _ ^)2^-i) 



4vr3/2i z3/2 

By Proposition 2.2 in |Ko2j . the latter expression is equal to exp(6*) A(;U, i/) 
whence (15.11). 



To prove (j5.3p and (j5.4p . we proceed similarly as in the last section. First of all, 
we rewrite the integrand as 

h[z)h[z) 

(1 - z)2 2{n+m)/2 ' 
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where now 



fiiz) :=exp (^-(4iV+ i(;x + z.)7Vi/3)i±|+rz^ • exp 



and 



/2(z) := exp i-w{z)) ■ la iw{z)) . 
Similarly as in the previous section, we find have 



=exp -m 



l-z 



+ b*z + 0{N-'^) 



l-z 



To prove ()5.3p . we use Taylor expansion. Straightforward calculations show that, 
for \u\ < a, 



h{l{u/N^/^)) 



f2{l{u/N^I^)) 



exp 



(_(l_i^)Ar2/3_i^ri/3_i+^(l_i^)3 



V^167riV4/3/(i 
(1 - 7(VA^i/3))2 ^ _ i^)2/^2/3 ^ o^(Ar-V3)^ ^ 

(7(n/iVi/3))('"+'^)/2 = exp ((1 - iu)N^^^ + iiV^/^ + i) (l + Oa(iV"^/^) 

For the second approximation, we have used the observation that for |m| < a, 
u;(i?e™/^) = m'^/^/il - iu){l + Oa{N-^l^)) , as well as the asymptotic 
approximation (13. 7p for the modified Bessel function. Putting it all together, 
the highest-order terms cancel out, and we find that the left-hand side of (jS.Sp 
is asymptotically given by 



exp(6*) exp (^(1 - iuf + - iu) - ^(// - z^)7(l " iu)) 



(1 - m)3/2 



du , 



which establishes (|5.3p . 

It remains to show ()5.4p . To this end, we will first show that 

MRe'') 



0{1). 



(5.5) 



Write y/Re^ = re''^, where r = y/l - l/iVV3^ g [-7r/2, -F7r/2]. Then, since 



Re 



1 - re*'^ 
1 + re*"^ 



1 - 



re 



.lip 



1 + + 2r cos 99 
1 



1 - fre*'^)2 



< 



iVV3, 
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the proof of (|5.5p is reduced to showing that 

expf-4A^i-^^'^ 



(_Rgjt)(m+n)/2 



Now, for sufficiently large N G'N, 



Re 



1 - re''^\ 
1 + re*"^ / 



1 



Ar-i/3 



1 + + 2r cos 99 2r(l + cos f) + {1 ~ r)^ ' 
where Cq denotes a constant depending only on a which need not be the same 



at each occurrence. It follows that 



exp ( -4iVi 



< 



exp 



< exp 



+ Ar2/3 , 1 ^1/3 , n] 

2r(l+cos</.) + (l-r)2 + 2^"^ + J 

' -47V2/3 + 2r (1 + cos ^)7V2/3 + ^(1 + cog y,)7vV3 + 
2r(l + COSV3) + (1 - r)2 



(Here some uniformly bounded terms have been absorbed into the constant Ca-) 
Since r = \/ 1 — N^^/^ < 1 — ^N^^^^, the numerator is bounded above by 

_4Ar2/3 + 2(1 + cos ip)N'^/''^ - (1 + cos ip)N^/^ + r(l + cos ip)N^/^ + Ca < Ca , 

and (jS.Sp is proved. 

Furthermore, the same arguments as those leading to (j4.8|) yield 



\f2iRe' 



(5.6) 



Hence, by (15. 5p . (15. 6p . and symmetry, in order to complete the proof of (|5.4p . 
it remains to show that for a > aQ{6), 



1 

2^ 



1 



dt < 5N^I^ . 



,a/m/3 |l-i?e^*|3/2 

But this can be proved in the same way as (j4.1ip . 

With essentially the same proof, we obtain the following result: 



(5.7) 

□ 



Theorem 5.2. Let f denote the second-order correlation function of a real sample 
covariance matrix satisfying our standing moment conditions. For any a E N, 
1/ G M, setting 

ZN{fi, v) ■= (im^ + A{n + u)N^/^ + ni^N^/^y^^ exp (^-AN - ^{fi + u)N^/^ 
we have 

f{n, m; AN + 2-^/^ fiN'^/^ , AN + 2^/3i/ArV3) \ 



hm 4 (2^/^,2^/31 

AT— >cxj \ 



n! ml 



exp(6*) A(/i, i^) , (5. 



where b* := {b — 3) and A is defined as in ()3.10p . 
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6. ASYMPTOTICS AT THE HARD EDGE OF THE SPECTRUM 

In order to zoom in at the hard edge of the spectrum, we have to make the 
replacements 

IJ. ;u/4iV and u u/AN , 
where G We then have the fohowing result: 

Theorem 6.1. Let f denote the second- order correlation function of a complex 
sample covariance matrix satisfying our standing moment conditions. For any 
Q G N, /U, 1/ G M+, setting Z^ilJ', v) ■= {fii' /N"^)"^"^ , we have 



lim 



= exp{b*)Sa{fJ',^) , (6.1) 



where b* := 2(6 — |) and Jq, is defined as in ([3 
Proof of Theorem \6.1[ By (|3.6|) , we have 

N-^ZNifJ.,!^)- 



r-1 r,_j.. .,^ /(n,m;^/iV,i//Ar) 



n! m\ 

1 



exp [-{ii + v)N-^ ■ + • lo^iwiz)) 

2TriN (1 - z)2 . z(n+™)/2 ~ ' ^^-^^ 

where now 

w(z) := 2^N^^ ■ . 

1 — z 

For the radius R we make the same choice as for the bulk of the spectrum, 
namely R = 1 — 1/N. Similarly as in the preceding sections, we reduce the proof 
to the following two claims: 

On the one hand, we will show that for any (fixed) a > 0, 

1 r exp (^-(/i + zy)iV-i • + b*z^ ■ Iaiw{z)) 

2^ iy (1 - zy ■ z("+-^y^ T ' "^^"^ ' ^^-^^ 

where 

1 /•+- exp ((1 - iu) - (f, + f)/(1 - ill)) I„{2^/(1 - iu)) 

On the other hand, we will show that for any d > 0, there exists some ao((5) > 
such that for a > ao{5), we have 



2-KiN J » (1 - z)2 • 2(n+™)/2 



< S (6.4) 



for all G N large enough. 

Once these two claims are established, the proof of Theorem 16.11 is completed 
using similar arguments as in the preceding section. ()6.3p and ()6.4p imply that 

1 



exp + v)N + b*z) ■ Ia{w{z)) 

> J(oo) , 



2TTiN J (1 - z)2 . ^(n+m)/2 
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where 
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J(oo) 



exp(6*) exp {z - (ii + u)/z) /„(2^/z) 



2m 



dz . 



1— ioo 



Substituting z = z/4:, dz = dz/4, and shifting the path back to the hne 1 + 
(which is easily justified by Cauchy's theorem), we obtain 



4exp(6*) exp {^z - 4(/i + i^)/z) Ia{8^/z) 



2m 



dz . 



7 

'1— ioo ^ 

Making the replacements /ii— >/i/4, i/h^i//4 and dividing by 4, we further obtain 
exp(6*) exp [\z - (/^ + v)lz) I^{2^/z) 



2m 



■ dz . 



1— ioo 



The subsequent Lemma [6.21 states that this is equal to exp(6*) JJa(^, '^)! thereby 
proving (16.ip . 



To prove (j6.3p . we proceed by Taylor expansion once more. A major difference 
to the previous situations is given by the fact that, for \u\ < a, 

u;(i?e»/^) = 2^N-^ • -^^^^—jjj = 2^ • (1 + Oa{l/N)) 



1 - Re^^l^ 
remains bounded. Since, for \u\ < a, 



1 — iu 



(1 - Re'-^/^f = (1 - iuf/N^ (1 + Oa{l/N)) , 

^^gWAf){n+m)/2 ^ exp(-(l - iu)) (1 + Oa{l/N)) 



we obtain 



1 I- exp(-(/i + zy)Af-i-^ + 6*z) 

z 



2mN Jy 
exp (6 



(1 - z)2 . 2{n+'n)/2 

f^" exp ( - + - in)) (2^/(1 - iu)) 
L (l-^.)VA^2-exp(-(l-.n)) i^ + o.il))du, 



27r7V2 
whence (|6.3p . 

(16. 4p follows immediately from the estimates 



< N. 



1-z 
and, for a > ao{5), 



1-z 



<N, \Iaiz)\< exp(Re z) 



^(n+m)/2 



1 



11 - i?e**P 



< 5iV. 



The last estimate is proved similarly to ()4.1ip . 



□ 
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Lemma 6.2. For any a G N, /i, G M+, 

1 exp {\z - (/^ + u)/z) 



27rz Ji^ioo z"^ 

Proof. We start from the following integral representation for the product of two 
Bessel functions (see e.g. p. 281 in [Erj): 

Differentiating with respect to x and y, we obtain 

1 /■^+^°° exp(-(x^ + y^)/z) {-2xU2xy/z) + 27/J4(2xj//z)) 



2vri ./i„joo ' 2^ 

= 2Vij'MxVt)Ja{2yVi) 



and 



1 /■^+^°° exp(-(x^ + y^)/z) i-2yU2xy/z) + 2xl'^{2xy/z)) 



2vr« Ji-ioo ■^^ 

= 2VtJa{2xVi)J'^{2yVt) 

respectively. Multiplying the former equation by x and the latter equation by y, 
taking the difference and dividing by 2(x^ — y^) yields 

J_ /•^+^°° exp(-(x2 + y^)/z)U2xy/z) 

_ Ja{2xVi)2yVt J^{2yVi) - 2xVt J^{2xVt) Ja{2yVi) 
~ 2(x2 - y2) • 

The assertion now follows by letting x := ^/JI, y := y^, and t := 1/4. □ 
The following result can be proved completely analogously to Theorem 16.11 



Theorem 6.3. Let f denote the second-order correlation function of a real sample 
covariance matrix satisfying our standing moment conditions. For any a G N, 
li^v ^ setting Zi^{ijl,v) := (hv/N'^)"^'^ , we have 



lim (l/16)iV-^Z7vW4,z//4) 



2^ / /. /.N f{n,m;fi/4N,u/AN) 



N^oo V n\m\ 

= exp(6*) JlQ,(/i,zy) , (6.5) 

where b* := (6 — 3) and Jq is defined as in ()3.1ip . 
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